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Abstract 

In this note we give a family of planar polynomial differential systems 
with a prescribed hyperbolic limit cycle. This family constitutes a corrected 
and wider version of an example given in the work 1 . The result given in 
this note may be used to construct models of Lienard differential equations 
exhibiting a desired limit cycle. 
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1 Introduction and statement of the main result 

Our purpose in this work is to give a family of planar polynomial differential systems 
of the form: 

x = P{x,y), y = Q{x,y), (1) 

for which an explicit expression of a limit cycle, that is, an isolated periodic orbit, 
can be given. We assume that P{x, y) and Q{x, y) belong to the ring of real polyno- 
mials in two variables y], and we will always assume that P{x,y) and Q{x,y) 
are coprime polynomials. We denote by d = maxjdegP, degQ} and we say that d 
is the degree of system ((T)). 

In the work [1^ a family of planar polynomial differential systems like is 
studied and the existence of an explicit limit cycle is pretended to be given. The 
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author of pQ gives a family of systems of the form with a prescribed invariant 
algebraic curve. This curve f{x,y) = has an oval surrounding the origin of 
coordinates. However, in there is no proof of the fact that the oval of f{x, y) = 
is an isolated periodic orbit, that is, a limit cycle. It is stated as obvious. We have 
been able to weaken the hypothesis appearing in jT] , getting a bigger family of planar 
polynomial systems, and we have been able to show that the oval of f{x,y) = is 
a hyperbolic limit cycle. 

Theorem 1 We consider a polynomial p{x) such that p{xe) = and p{xd) = for 
some Xe < and Xd > 0, p'{xe) 7^ and p'{xd) 7^ 0. We assume that p{x) > 
for all X in the interval {xe,Xd)- We consider another polynomial q{x) satisfying 
p{x)q{xY 7^ 1 for all x G {Xe,Xd) and q\x) 7^ for all x G {Xe,Xd)- Then, the 
algebraic curve given by f{x,y) = with f{x,y) := {y — p{x)q{x)Y — p{x) has an 
oval in the band x^ < x < Xd which is a hyperbolic limit cycle for the following 
system: 

y, 

(^q{x)p'{x) +p{x)q'{x)^ y - {p{x)q{xY - l) . 



X = 

y = 



We note that a system of the form Q can also be viewed as an autonomous 
Lienard differential equation: 

X + f{x)x + g{x) = 0, (3) 

where f{x) and g{x) are the polynomials given by: f{x) := —3q{x)p'{x)/2 — 
p{x)q'{x) and g{x) := {p\x)/2) {p{x)q{xY — 1). Therefore, Theorem [T] may be used 
to construct models of Lienard differential equations exhibiting a hyperbolic limit 
cycle. 

One of the most famous Lienard differential equation is called van der Pol equa- 
tions and it appears when studying the vacuum-tube circuits. This particular equa- 
tion Q has f{x) = ^{x"^ — 1) and g{x) = 1, with /i G M, and it exhibits a unique 
hyperbolic limit cycle surrounding the origin. This limit cycle is shown to be non- 
algebraic in the work of Odani 0. We are not considering van der Pol's equation 
since the systems described in Theorem ^ always exhibit an algebraic limit cycle. 
The systems given in (j2)) are examples of Lienard equations with an algebraic and 
hyperbolic limit cycle. 

The same kind equations are studied in the work P] , but under other hypothesis 
for the polynomials p{x) and q{x), and the author of pretends to state the exis- 
tence of a limit cycle. The conditions for the polynomials p{x) and q{x) appearing 
in pp are: p{x) is an even polynomial, p{0) > 0, there exists a value X > such that 
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p{X) = 0, p{x) < for all x > X, q{x) is an odd polynomial and g(0) = 0. All these 
conditions are contained in the ones that we assume in Theorem ^ However, the 
authors noticed that in the work Pl the condition p{x)q{x)'^ ^ 1 for x G (— X, X) 
does not appear and it is not implied by the other hypothesis. As we will see in the 
proof of Theorem HJ condition p(x)g(a;)^ 7^ 1 for a; G (xe,^^) is necessary to have a 
limit cycle and it cannot be avoided. 

We remark that each system of the family Q has an algebraic limit cycle, which 
is an oval of the real algebraic curve f{x,y) = 0, and it may have other limit cycles 
which are not taken into consideration. These other limit cycles can be contained 
in f{x, y) = OT not. If they are contained in an invariant curve, we can treat them 
with the same methods described in this note. For instance, it can be shown that 
the system Q with p{x) = (1 — a;^)(4 — x^)(9 — x^) and q{x) = x/100 has 3 hy- 
perbolic limit cycles all of them contained in the corresponding invariant algebraic 
curve f{x, y) = 0. 

In order to prove Theorem ^ we need some results relating ovals of curves of the 
form f{x,y) = with the fact that the oval of such a curve is a limit cycle. These 
preliminary results are stated in Section |21 Once we have stated these previous 
results, we prove Theorem [T] in Sectional 



2 Preliminary results 

We are considering limit cycles which are contained in a real curve f{x,y) = 0, 
which does not need to be algebraic. This fact leads us to the definition of invariant 
of a system (P). 

Definition 2 Let us consider an open setU C and aC^iU) real function denoted 
by f{x, : W C — >• M. We say that f{x, y) is an invariant for a system if 

Of df 
P{x,y) -^{x,y) + Q{x,y) —{x,y) = k{x,y) f{x,y), (4) 

with k{x,y) a polynomial of degree lower or equal than d — 1, where d is the degree 
of the system. This polynomial k{x,y) is called the cofactor of f{x,y). 

In case that f{x,y) is a polynomial we say that f{x,y) = is an invariant al- 
gebraic curve for system (HJ. We notice that if f{x,y) is an invariant of sys- 
tem and f{x,y) = defines a curve in the real plane, then the function 
P{x,y) (df/dx) + Q{x,y) {df /dy) equals zero on the points such that f{x,y) = 0. 
This fact implies that the real curve f{x,y) = is formed by orbits of system (0). 
In particular if f{x,y) = contains an oval without any singular point of system 
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(P), this oval is a periodic orbit of system ([T]). 



As well as invariant curves, the other objects taken into consideration in this 
paper are limit cycles. A limit cycle of system is an isolated periodic orbit. Let 7 
be a limit cycle for system (P). We say that 7 is stable if there exists a neighborhood 
such that all the orbits starting in it have 7 as cj-limit set. We say that 7 is unstable 
if there is a neighborhood such that all the orbits starting in it have 7 as a-limit 
set. There might be limit cycles which are neither stable nor unstable. These limit 
cycles have a neighborhood such that in the interior of the limit cycle all the orbits 
have 7 as tu-limit set and in the exterior of 7 all the orbits have 7 as a-limit set. Or 
the other way round: the orbits of the interior have 7 as a-limit set and the orbits 
in the exterior have 7 as tu-limit set. In this case, we say that 7 is semi-stable. Any 
limit cycle 7 of a system (H)) is either stable, unstable or semi-stable as it is stated 
in m. 

A classical known result, given in the book of Perko let us distinguish the hy- 
perbolicity of a limit cycle. If we consider 7(t) a periodic orbit of system of period 
T, we may compute the finite value given by the following integral Jg^div(7(t)) dt, 
where div(x, y) = {dP / dx)-\-{dQ / dy) is called the divergence of system (H)). It can be 
shown that if div(7(t)) dt < 0, then 7 is a stable limit cycle, if div(7(t)) dt > 0, 

then 7 is a unstable limit cycle and if div{'y{t)) dt = 0, then 7 may be a stable, 
unstable or semi-stable limit cycle or it may belong to a continuous band of cycles. 
When the quantity div(7(t)) dt is different from zero, we say that the limit cycle 

7 is hyperbolic. We notice that if dw{'y{t)) dt 7^ 0, then the periodic orbit 7 is a 
limit cycle (either stable or unstable). We are going to use this property to ensure 
that a periodic orbit is a limit cycle, that is, that it does not belong to a continuous 
band of cycles. 

We relate limit cycles with invariants in the following way. We assume that we 
have a periodic orbit 7 of system which is given in an implicit way, that is, there 
exists an invariant curve f{x,y) = such that 7 C {{x,y) \ f{x,y) = 0}. In order 
to have a smooth curve f{x,y) = defining the periodic orbit, we will assume that 
V f{p) 7^ for any p G 7, that is, the gradient vector of /(x, y) is different from zero 
in all the points of 7. Then we have the following result stated and proved in [2]. 

Theorem 3 Let us consider a system ^ and'~f{t) a periodic orbit of period T > 0. 
Assume that f : U M."^ —>■ is an invariant curve with 7 C {{x,y) \ f{x,y) = 0} 
and let k{x,y) be the cofactor of f{x,y) as given in We assume that V/(p) 7^ 
for any p G 7. Then, 
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Hence, we have an alternative way to compute the value div{'-f{t)) dt. In the 
family of planar polynomial differential systems which we are considering, that 
is, the one described in Theorem Q we will not be able to directly compute the 
value Jq div {•y (t)) dt. This is due to the fact that we are not considering a fixed 
system with a concrete periodic orbit, but a family of systems each one with a 
different periodic orbit and, thus, the expression of the integrand is too general to 
be manipulated. By Theorem |2l we can also compute the value k{'j{t)) dt but 
this integral is as much difficult as the previous one. That's why we are going to 
use the fact that for any w G M: 

div {-f (t)) dt = div {-f (t)) dt + w i div {-f (t)) dt - k{j{t)) dt 
Jo \Jo Jo 

The integrand in the right hand side of this equality will be chosen strictly positive 
or negative in all the interval of integration for a suitable value of w. Therefore, the 
value of the integral will be different from zero. Using these steps, we will be able 
to prove that the oval of the invariant curve described in Theorem ^ is a limit cycle 
of the corresponding system. 



3 Proof of Theorem [T] 

In order to prove this theorem, we first show that f{x,y) = 0, where f{x,y) : = 
{y —p{x)q{x))'^ —p{x), is an invariant algebraic curve of system Q. Straightforward 
computations show that: 



^1 + 

dx 



^ q{x)p'{x) + p{x)q'{xyj y - ^^-^ {p{x)q{xf - l) 



dl 
dy 



q{x)p'{x)f{x,y), 



and, thus, we have that f{x,y) = is an invariant algebraic curve for system Q 
with cofactor k{x,y) := q{x)p'{x). 

Since p{xe) = p{xd) = for the values Xg < and > and p{x) > for 
X G {xe,Xd), we deduce that f{x,y) = has an oval in the band Xe < x < Xd 
surrounding the origin of coordinates, which can be parameterized in two parts by: 



x{t) = T, y±{T) = p{T)q{T) ± Vp(r), (6) 

with r G {xe.Xd)- We are going to prove that this oval does not contain any 
singular point of system and then, we will have that it defines a periodic orbit 
of the system. The singular points of system Q have coordinates of the form (a, 0) 
where the value a is a root of the polynomial p'(x)(p(a;)g(x)^ — 1). We have that 
/(a,0) = p{o){p[a)q{a)'^ — 1) and {p{a)q{aY — 1) is different from zero in all the 
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closed interval a G [xejX^] by the hypothesis that 7^ 1 for all x G {xe,Xd) 

and p{xe) = p{xd) = 0. 

Here we notice that the assumption p{x)q{x)'^ 7^ 1 for x G {xe,Xd) is necessary 
for the oval of f{x, = to be a limit cycle. In yy, this assumption is not given. 
We notice that an oval of an invariant algebraic curve of a system may contain 
singular points of the system, and in such a case, it is not even a periodic orbit. 

Since x^ and Xd are simple zeroes of p{x) and p{x) > for all x in the interval 
{xe,Xd) we have that there is no singular point of system Q on the oval given by 
f{x,y) = and parameterized by ©. From this fact and that f{x,y) = is an 
invariant algebraic curve of the system, we deduce that this oval is a periodic orbit 
of system (jT)). We denote this periodic orbit by 7 for the rest of the proof. We 
note that we do not know the parameterization of 7 as explicit solution of system 
Q, that is, we do not know the periodic function 7(t) := (7i(t), 72(i)) such that 
d'yi(t)/dt = 72 (t) and 

"^^'^^^ - ^??(7iWy(7iW)+p(7i(t))9'(7i(t)))72(t) 

{p{ii{tMii{t)r - 1) , 



dt \2 

P'(7iW) 
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for all t G M. We do neither know its period T > but we have been able to show 
its existence by using the invariant algebraic curve f{x,y) = and its properties in 
relation with system 

Finally, we need to prove that the periodic orbit 7 is a hp^erbolic limit cycle. To 
do so, we are going to show that the value of the integral div(7(t)) dt is different 
from zero. Since we do not know 7(t) nor the period T, we use the parameteri- 
zation of the oval 7 given in In order to get the correct sign of the integral 

div(7(t)) dt, we need to know the sense of the flow over 7. We take the point of 
coordinates (xrf,0), which belongs to 7, and we have that the vector field defined 
by system ^ on that point is {0,p'{xd)/2) because p{xd) = 0. Since p{xd) = 0, 
p'{xd) 7^ and p{x) > in the interval {xe,Xd), we deduce that p'{xd) < 0. Hence, 
the sense of the flow over 7 is clockwise. We can write the following equality, using 
the parameterization 



div(7(t))c^t = / -div(x(r) y,(r))^^^ div(x(r) 



^div(x(r),?/+(r)) div(a;(r), ?/_(r))\ 



T 



We note that the divergence of system (j2|) is div(x,?/) = 3g(x)p'(x)/2 + j9(a;)g'(x), 
and substituting this expression in the former equality, we get: 

div 7 t rft = / — '^^ ' Ll_ ' dr. 
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This integral is well defined because we are assuming that {p{r)q{ry — 1)p{t) is 
different from zero for r G {xe,Xd)- However, we are not able to distinguish if its 
value is positive, negative or zero. By using the same reasonings, we can write the 
following equality: 

r^MiMi) 

A, (p(tMt)2-1)VpM 
Using Theorem El we have that, for any value of w G M: 

div(7(t))cit = / div(7(t))dt + w ( / div {-f (t)) dt - j k{-f{t))dt 

dr. 



\Jo Jo 

-{w + 3)g(r)p'(r) - 2(1 + «;)p(r)g'(r) 



(p(r)g(r)2-l)VR^ 



Taking w = — 3, we get: 



The hypothesis of Theorem^ on the sign of the polynomials p{x) and q{x) in the 
interval x G {xe,Xd) are p{x) > 0, p{x)q{x)'^ ^ 1 and q'{x) ^ 0. Therefore the 
integrand of the right hand side of ((7j) is strictly positive or negative in all the 
interval r G {xe,Xd)- We deduce that the value of the integral cannot be zero and, 
hence, the periodic orbit 7 is a hyperbolic limit cycle as we wanted to show. 

We would also like to characterize if this limit cycle is stable or unstable, so we 
are going to study the sign of the integrand in the right hand side of (|7j). Since 
p{x)q{xY 7^ 1 for X G (xg, Xd) and p{xd)q{xd)'^ — 1 = — 1 (because p{xd) = 0), we de- 
duce that p{x)q{xy — 1 < for x G {xe,Xd)- Therefore, we have that the integrand 
in the right hand side of (|7j) is strictly positive if ^'(r) < for all r G {xe,Xd) and 
strictly negative if q'i^r) > for all r G {xe,Xd)- We can state that the hyperbohc 
limit cycle 7 is stable if g'(0) > and unstable if g'(0) < 0. | 

We also note that in the work pp , the expression of an example of a more general 
limit cycle for a family of planar polynomial differential systems is pretended to be 
given. In fact, we are going to show that in the case that the oval of this example 
is a limit cycle, we are in the same family of systems as written in that is, the 
one described in Theorem ^ 

In PP the following planar polynomial differential system is given as an example 
of a more general family of systems with an explicit limit cycle. 

X = y, 

y = {p'{x) [{m + r)h{x)p{xy~^ + (m + l)q{x)~\ + h'{x)p{xY ^g") 
+p{x)q'{x)} y — mp{x)p'{x) ([h{x)p{xy'~^ + <l{x)Y — p(a;)^"'~^j , 
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where p{x) is an even polynomial, q{x) and h{x) are odd polynomials, m = n + 1/2, 
n is an integer number with n > and r is an integer number with r > 2. Moreover, 
it is assumed that p{0) > and there exists a value X > such that p{X) = 0, 
p{x) < for all X > X and g(0) = 0. 

Some straightforward computations show that system (jHJ exhibits the invariant 
algebraic curve /(x,?/) = with /(x, := {y — h{x)p{xY — q{x)p{x)Y —pi^xY"^^^ and 
with cofactor k{x,y) := {2n + l)p'{x) [h{x)p{xY~^ + q{x)]. We have that f{x,y) = 
has an oval in the band —X < x < X which can be parameterized by: 

x{t) = T, y{T) = /i(r)p(r)" +p(r)g(r) ±p(r)"V^. 

In order to show that this oval is a periodic orbit of system (jH)), we only need to 
show that it does not contain any singular point of the system. The singular points 
of system (jH)) in the band |x| < X are of the form (a, 0) where a is a root of the 

polynomial p{x)p'{x) (^[h{x)p{xY^^ + q^x)]^ — j , because m = n + 1/2. We 

notice that, unless n = 0, the points with coordinates {—X, 0) and {X, 0) are singular 
points of the system which are contained in the oval of f{x,y) = 0. Therefore, if 
n > 0, we have that the oval of f{x, y) = cannot be a limit cycle. If n = 0, we can 
consider the polynomial q{x) := q{x) + h{x)p{xY~^ and we have that system (jH)) 
coincides with system ^ with polynomials p{x) and q{x). Therefore, this is not an 
example of a more general limit cycle. 

References 

[1] M.A. Abdelkader, Relaxation oscillators with exact limit cycles. J. Math. 
Anal. Appl. 218 (1998) 308-312. 

[2] H. GlACOMlNl AND M. Grau, On the stability of limit cycles for planar 
differential systems, J. Differential Equations 213 (2005), 368-388. 

[3] K. Odani, The limit cycle of the van der Pol equation is not algebraic, J. 
Differential Equations 115 (1995), 146-152. 

[4] L. Perko, Differential equations and dynamical systems. Third edition. Texts 
in Applied Mathematics, 7. Springer- Verlag, New York, 2001. 



8 



